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We investigate the lattice index theorem and the localization of the zero-modes for thick classical 
center vortices. For non-orientable spherical vortices, the index of the overlap Dirac operator differs 
from the topological charge although the traces of the plaquettes deviate only by a maximum of 
1.5% from trivial plaquettes. This may be related to the fact that even in Landau gauge some links 
of these configuration are close to the non-trivial center elements. 
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I. INTRODUCTION 



At large separations static color sources of non-vanishing N-ality give rise to gluonic flux tubes and a linearly 
rising potential. This behavior is caused by the non-perturbative QCD-vacuum which compresses the gluonic flux. 
d ^ Center vortices P, 0, H 0, H, @] , quantized magnetic flux lines, play an important role for the confinement of color 
I ' ' charges, as has been shown by numerical calculations 0,0]. In addition, simulations have indicated that vortices could 
C^j also account for phenomena related to chiral symmetry, such as topological charge and spontaneous chiral symmetry 
^ breaking (SCSB) |, 0, [HI, [ll • 

1^ I ■ It is well-known that these non-perturbative features of QCD are intimately linked to the properties of the low- lying 
spectrum of the Dirac operator. The topological charge of a gauge field equals the index of the Dirac operator, while 
' the chiral condensate, the order parameter for SCSB, is proportional to the spectral density of the near-zero modes. 
With the advent of the overlap operator [l^ [3] , the fundamental problems of investigating chiral features on the 
lattice have been overcome. The overlap operator has an exact chiral symmetry 15], implements a lattice version of 
\ the index theorem [H] , and may even be used for the definition of a local topological charge density [3| . 
• In this paper, we report on our calculations with the overlap operator applied to thick classical center vortices. We 
investigate the localization of zero-modes with respect to the position of the thick vortices and find an interesting 
discrepancy in the topological charge determined by different methods. 



II. STRATEGY 

We work with SU(2) gauge fields on lattices with Ng sites in the x-, y- and z-direction, and Nt sites in the t-direction. 
In four dimensions vortices form closed two-dimensional surfaces. We investigate thick classical center vortices in the 
shapes of planes (closed by lattice periodicity) and spheres. The details of the individual vortex types will be discussed 
along with our results. 

We compare different definitions of the lattice topological charge: 

(1) The topological charge determined from the shape of P-vortices, that means from vortex intersections and 
writhing points according to ref. Q . 

(2) The index of the overlap Dirac operator [l^. 

According to the Atiyah-Singer index theorem, which holds for continuous fields, the topological charge is given 
by the index 

indD[A] = - n+ = Q (1) 
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where n_ and n+ are the number of left- and right-handed zero modes [H, [ill • It has been shown that 
this relation is also valid on the lattice for the overlap Dirac operator: its index coincides with the discretized 
gluonic definition in the continuum limit [l9l |. 

The overlap Dirac operator is defined by [l^l 



D = Ti[l + l^e{H+)] (2) 



Here, e is the sign function, 

H+ =-i^D.^{-mQ), (3) 
and Dw is the usual lattice Wilson Dirac operator with mass — toq (we use mg — 1.0). 

The overlap operator can be used to define the topological charge only on so-called "admissible" gauge fields. 
This restriction assures that iJ^ has no zero eigenvalues so that the sign-function is well-defined. A sufficient, 
but not necessary condition for this is [U, 

tr(l - U^,) < 0.03 (4) 

The diagonalization of the overlap operator yields the corresponding Dirac zero modes with well-defined chirality 
properties. We further look at the distribution of the scalar fermionic density 

c,d 

where the summation indices c and d refer to color and Dirac indices of the eigenvectors v, respectively, and the 
left- or right-handed chiral densities 



P±{x) = ^(^Tcd— ^(^)crf'> (6) 

'd.d 



d.d' 



where only zero-modes of the corresponding chirality contribute. When the zero modes of a specific chirality 
are degenerate, it is useful to investigate the square root of the sum of the scalar densities of these modes. This 
average density respects eventual symmetries of the field configuration. 

(3) The integral (sum, on the lattice) of the gluonic charge density q{x) — tr(jFp^J^p^) in the "plaquette" 
and/or "hypercube" definitions on the lattice, see ref. [Ij, [23|. Since Monte-Carlo configurations are in general 
too coarse, these definitions are usually only applied after cooling. 

In the continuum, since the topological charge density q{x) is the total derivative of the topological current fc^ 



q{x) = 5pfcp, 

2 
^3 



(7) 



one can apply Gauss' theorem to transform the expression for Q into a surface integral, 

Q = J '^^^ '^^^^ ^ j' ^'^m'^^M' (8) 

provided the gauge field is smooth. Assuming a sufficiently fast decaying field strength for a;^ — > cx), the gauge 
field at infinity can be written as a pure gauge 

^^(a:2-^c5o) = i9^f7f]t. (9) 

It can be shown that the surface integral Eq. ([8|), and hence Q for smooth fields, measures the winding number 
of ri, which maps the group space SU(2) onto the space-time surface S"^ at infinity. 
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Geometry Topological charge density Fermionic density 




FIG. 1: Plane vortices on a 12^-lattice in xy- and zt-plane intersect in four points giving topological charge and fermionic 
density. 

III. DESCRIPTION OF CONFIGURATIONS AND RESULTS 
A. Plane vortices 

The construction of planar vortices is explained in |2Q]. We use periodic boundary conditions where plane center 
vortices always come in pairs. All three definitions of Q yield identical results for all configurations containing only 
plane vortices. 

As an example, we show in Fig. [T]some diagrams for two orthogonal pairs of plane vortices on a 12^* lattice, which 
intersect in 4 places. The left picture shows the position of the P-plaquettes after going to maximal center gauge 
and center projection. One vortex pair extends in x-y-planes, the other pair in z-t-planes. The short pieces of lines 
attached to the vertical line symbolize the extent of the z-t-vortices in the time-direction. The diagram in the middle 
of Fig. [T] shows an equi-density surface of the topological charge density determined in the plaquette definition and 
the right diagram shows an equi-density surface of the scalar fermionic density. 

Each of the four intersection points gives rise to a lump of topological charge Q = ±^@. The sign of the contribution 
at a given intersection point can be changed by a flip of the orientation of the vortex surface, this means by a transition 
from links U^{x) creating one of the vortices to U'^{x) in the region of the intersection. After abelian projection a 
smooth transition of the vortex orientation in color space produces a monopole line. In this way one can produce 
configurations with two vortex pairs of topological charge Q = 0, ±1 or ±2. In all these cases the lattice index theorem 
is fulfilled and our results agree with the analytical solution for the zero modes presented in p7j . 



B. A non-orientable spherical vortex 

A non-orientable spherical vortex can be constructed with an appropriate "profile" -function a{\r — ro|) by time-like 
links which form a hedgehog in one time-slice of the lattice, where fo is the center of a spatial sphere 

I 1 elsewhere. 

A characteristic property of a hedgehog configuration is the agreement between the color direction n and the spatial 
direction of r — ro 

n(r,t) = ^^, (11) 
\r - ro\ 

see left diagram in Fig. [2l For the choice (flO)) of the links the time-like Wilson lines (holonomies) agree with the 
time-like links in the t = 1 time-slice. 

At r = fo and at the largest possible distances from rg on a periodic lattice the direction of the color-vectors 
pT|) are undefined. This does not lead to a singularity of the gauge field if the profile-function a(r) is appropriately 
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chosen such that the Wilson hnes at these special points become center elements of the gauge group. A good choice 
is a(r) = a-{r) with 

r < i? - f 

R-^<r<R+^ (12) 
i? + f < r. 

As shown in the right diagram of Fig. [2] for a 40"^ x 2 lattice the phase a of the links in the first time slice varies from 
TT at the center of the sphere to outside of the vortex core of thickness A, chosen as A = 18. 

The check that this configuration is a vortex is done with maximal center gauge and center projection and results 
in a P-vortex forming a lattice representation of a sphere of radius 10 in the first time slice. The color structure of 
the thick vortex which is symbolically indicated in the left diagram of Fig. [3] leads to a monopole loop on a great 
circle of the P-vortex after maximal abelian gauge and abelian projection. The direction of the loop depends on the 
U(l) subgroup chosen as abelian degrees of freedom. For the subgroup defined by the Pauli matrices cti, a2 or ag, the 
monopole loops are in the y-z-, z-x- and x-y-plane, respectively. This is indicated schematically in Fig. [3] with three 
colors. 






FIG. 3: The color structure of the vortex surface (1) leads to monopole lines after abelian projection (r). 



The P-vortex surface which is a closed surface in dual space consists only of the duals of space-time plaquettes. 
There are no vortex intersections and no writhing points. Therefore, the topological charge Q determined from the 
P-vortex surface vanishes, Q = 0. The traces of all plaquettes are close to unity, tr(l — Uf^,^) < 1 — cos = 0.015, 
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FIG. 4: Phases of Wilson lines along a line through the center of a spherical vortex (1), normalized scalar densities for all 
zero-modes of one chirality in the planes with z = 20 close to the center of the sphere (r). The density vanishes in the region 
where the phase of the Polyakov loop is tt. Both chiralities and both time-slices on the 40^ x 2-lattice give the same picture. 

compare the link-phases in Fig. [2l The topological charge in the plaquette definition or in the hypercube definition can 
be used safely without cooling. According to Eq. ([T0|) the only non-trivial links are time-like, therefore all space-space 
plaquettes are trivial, which confirms the vanishing of the topological charge. 

This configuration has three zero-modes [2§| of positive chirality and four zero-modes of negative chirality. Nor- 
malized densities of positive and negative chirality modes have the same spatial distribution. 

For the spherical vortex with center at (19.5,19.5,19.5), we show the scalar density in the plane z = 20 in Fig. 21 The 
non-trivial links connect z = 20 and z — 21, and the densities in both time-slices are identical. The zero- modes are 
localized mainly in the region of trivial Wilson lines (for |r — ro| ~ 20) and seem compressed in the regions where the 
mirror pictures of the spherical vortex on the periodic lattice approach each other. In the region where the Wilson 
lines arc non-trivial center elements (for |r — ro| around 0), the density vanishes. According to the index theorem the 
topological charge is given by 

Q = n_-n+ = l. (13) 

This value is in disagreement with the values determined above (Q = 0). But it agrees with the gluonic result which 
one gets after cooling. The cooling history is shown in Fig. [5l The left scale shows the action S — X]n(l ~ ^Tr[/n) 
in units of the one-instanton action S'inst and the right scale the value of the topological charge in the plaquette 
definition. The spherical vortex is not a minimum of the action, i.e. not an instanton. During cooling its action is 
first rapidly decreasing, whereas its topological charge starts with zero and reaches a value close to one after a few 
cooling steps. On the 40^ x 2-lattice the action does not really approach a plateau value. Since a plateau appears for 
a 12^-lattice this seems to be related to the strong asymmetry of the lattice. Hence one would expect to see a clear 
plateau on a symmetric, but finer 40* lattice, as well. The number of zero-modes is the same for various lattice sizes 
we used: 12'', and 12^ x 2 up to 40'^ x 2, increasing the spacial extent Lg in steps of 4. We conclude that cooling and 
overlap fermions give a different topological charge than the discretized integral J d'*a: tr(jF^yjFp^) even for very 
smooth field configurations. 

The hedgehog configuration is characterized by links in the center of the spatial sphere which are non-trivial center 
elements and Wilson lines L{x), x = {x, y, z), with a topological winding number = 1, 

^({i}) = ^ / d^re.,kTr{V,V,Vk), V, = . (14) 

The topological charge that we get in the cooling process and from the index theorem agrees with this winding number. 
In the continuum limit, the gauge field at the center of the hedgehog becomes singular. This singularity is due to the 
topological obstruction and can't be removed by a gauge transformation. With Landau-gauge, which we did for 100 
gauge copies on a 40^ x 2-lattice, it is not possible to get rid of links close to non-trivial center elements. For the 
lattice spacing a approaching zero, the gauge field for these links, Ai){x) ~ 'Kay,/{ga), diverges. In all gauge copies we 
found the same links with negative traces. The number of such negative links and the values of their traces are listed 
in Table [J 











10 20 30 40 50 60 70 80 90 100 
number of cooling steps 

FIG. 5: By cooling the topological charge rises from zero close to one for q = a_ (right scale) while the action S (in units of 
the one-instanton action Sinst) decreases slowly (left scale). The data is for a 40'^ x 2-lattice. 



number of links 


8 24 24 32 48 24 48 




-0.93 -0.68 -0.46 -0.30 -0.18 -0.08 -0.01 



TABLE I: Number of links with negative traces and trace values after Landau gauge for the non-orientable spherical vortex of 
Eq. dTO]) on a 40^ x 2-lattice, e.g. there are eight links with |tr(7^(x) = -0.93. 



C. Other configurations with spherical vortices 

An orientable spherical vortex we construct with the prescription 



exp {ia{\r~ fo|) |n,;(a;)| • ctJ i = 1, ^ = 4, 
1 elsewhere. 



(15) 



Here, 



ir-ro)i 

\r—ro\ 



is the absolute value of the coordinates of the direction vector ft. 



These links cover only 

half of the SU(2) group space and this half they cover twice. The time-like Wilson lines have topological winding 
number zero. This configuration has no zero-modes and all ways to determine the topological charge lead to the same 
result, Q — 0. With maximal Abelian gauge no monopole lines appear. With Landau-gauge, which we did for 100 
gauge copies on a 40^^ x 2-lattice, one can get rid of links close to non-trivial center elements. The number of such 
negative links and the values of their traces are listed in Table [III 



number of links 


8 24 24 24 


^tr(7p(x) 


-0.26 -0.15 -0.06 -0.008 



TABLE XL Number of links with negative traces and trace values after Landau gauge for the orientable spherical vortex of 
Eq. JlSl) on a 40^ x 2-lattice. 

A second type of a non-orientable spherical vortex we construct with the prescription (jlOp and the profile function 
a{r) = a+(r), 

R~^<r<R+^ (16) 
R+%<r 

The phases of the links in the time slice agree again with the phases of the Wilson lines in time direction and are 
shown in Fig. [51 In the center of the spherical vortex the temporal Wilson lines are trivial and at "infinity" they are 
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non-trivial center elements. This configuration gives only one zero-mode of negative chirality and none of positive 
chirality. The zero-mode is localized in the center of the spatial sphere and in both time-slices occupied by the vortex. 
Again, the scalar density prefers regions of trivial Wilson lines over those of non-trivial values. According to the index 
theorem the topological charge is -1 and agrees with the winding number (|14p of the hedgehog of Wilson lines. The 
topological charge computed from plaquettes is Q = while the plaquette traces differ by less than 1.5 % from trivial 
ones. Since this configuration is related to the configuration in Fig. |4] by a non-periodic gauge transformation, it is 
not surprising that again the different determinations of the topological charge do not agree. 




FIG. 6: Phase-diagram of the Polyakov loop (1) along a line through the center of a spherical vortex with the profile function 
a — a+ (1) and scalar density of the only fermionic zero mode (r) which has positive chirality. 

For two equal vortices of this type, a — «+, one in time slice t = 1 and another at < = 2 of the 40"^ x 2-lattice we 
get three zero-modes of positive chirality with a scalar density like the diagram in Fig. |4] and five of negative chirality 
with the densities shown in Fig. [T] 

Since the modes are degenerate, they can mix and the plots display the sum of the densities of all zero-modes. The 
left diagram displays the t — 1 slice, where only four of these modes contribute significantly, while the fifth mode 
leads only to the small peak in the middle. Conversely, the density in the t = 2 slice (shown in the right picture) is 
dominated by the fifth mode in the center, the other modes being of a magnitude invisible on this scale. 

The phase of the Polyakov line for this configuration is zero in the center of the spherical vortices and due to the 
27r-periodicity in the region between the mirror pictures on the periodic lattice. Again the zero-modes are localized 
in the region with trivial time-like Wilson lines. 




t=l t=2 

FIG. 7: Scalar densities of fermionic zero modes of negative chirality for two equal vortices with the profile function a = a+, 
in the time slices t = 1 and f = 2 of the 40"^ x 2-lattice. 
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A configuration with two vortices of different type in the two time-shces, one with a = a- at t = 1 and another 
with a = a_(_ at i = 2 does not show any zero-mode. In this case there is no discrepancy in the determinations of the 
topological charge. The phase of the Wilson line, which is the sum of the phases from the two vortices, is constant 
and equal to tt. 

A summary of the the number of zero-modes for the configurations described above with spherical vortices is given 
in Table ED 







profile 


a 






n+ 


n_ 






Q_ at t 


= 1 






3 


4 






a+ at t 


= 1 






1 





a_ 


at t 


= 1 and 


a+ 


a.t t — 


2 


3 


5 


«+ 


at t 


= 1 and 


Q_ 


a,t t — 


2 


5 


3 


a+ 


at t 


— 1 and 


a+ 


at t = 


2 









TABLE III: Number of zero-modes for various smooth field configurations with spherical vortices. The profile functions a± are 
defined in Eq. dHJ and (fTU)) . 

The above table can be summarized in the following empirical rule: Spherical vortices in slices (3D volumes) of the 
lattice contribute to the Dirac operator index with an integer given by the winding number (|14p of the corresponding 
Wilson lines, mapping the 3D volume of the slice to the SU(2) manifold of the Wilson lines. 

IV. CONCLUSION 

We have investigated various configurations with plaquettes close to unity. For configurations with a non-zero 
winding number of the Wilson lines we get disagreement between the topological charge determined by the plaquette 
and hypercube definition, and from the vortex intersections and writhing points on one side, and the topological 
charge from the index theorem and cooling on the other side. The transition in Eq. ([5]), from the 4D integral over 
the topological charge density, to the 3D integral of the SU(2)-valued map il, is only possible for smooth gauge fields. 
Vortex configurations with non-trivial winding numbers have singularities in the gauge fields which can't be removed by 
gauge transformations or smooth deformations of the field configuration. According to the vortex model of confinement 
vortices form closed two-dimensional surfaces which in the confined phase penetrate the whole space-time. In ref. (28j 
it was shown that P-vortex surfaces are mostly non-orientable. It is very likely that this non-orientability is related to 
a non-zero winding number of the gauge field and the artificial configurations which we discussed in this article reflect 
properties of those field configuration which saturate the path integral. The discrepancies which we have indicated are 
not due to the coarseness of the lattice, the traces of all plaquettes are close to unity, tr(l — U^^) < 1 — cos ^ = 0.015 
but probably due to singularities in the gauge field related to non-trivial mappings of Wilson lines. 
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